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We present theoretical study of a single disclination motion in a thin free standing hexatic (or 
smectic-C) film, driven by a large-scale inhomogeneity in the bond (or director) angle. Back-flow 
effects and own dynamics of the bond angle are included. We find the angle and hydrodynamic 
velocity distributions around the moving disclination, what allows us to relate the disclination 
velocity to the angle gradient far from the disclination. Different cases are examined depending on 
the ratio of the rotational and shear viscosity coefficients. 

PACS numbers: 05.20, 82.65, 68.10, 82.70.-y 
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CZ3 ' Introduction 
-f-i 

Physics of thin liquid-crystalline films has been a recurrent hot topic during the past decade because of their 
intriguing physical properties and wide applications in display devices, sensors, and for many other purposes. Hexatic 
1-^ I and smectic-C liquid crystalline films, we consider, belong to two-dimensional systems with spontaneously broken 
<—j . continuous symmetry. Therefore an essential role in the behavior of the films is played by vortex-like excitations 
(disclinations) . Defects are almost necessarily present in liquid crystals, and their dynamics plays a crucial role in 
the overall pattern organization. Early studies of defects focused on classifying the static properties of the defects 
and their interactions [1, 2]. More recently the focus has moved to examining the dynamics of defects (see, e.g., 
7-H ! [3] and references herein). Note that though in most practical applications of liquid crystals, such as traditional 

■ display devices, defects destroy an optical adjustment, and therefore are not desirable, there are novel display designs 
(bistable, multidomain liquid crystalline structures) just exploiting defect properties. 

Unfortunately, theoretical researches of dynamic characteristics of the films are in a rather primitive stage, in 
spite of the fact that experimental dynamic studies are likely to be more fruitful. It is largely accounted for by a 
complexity of dynamic phenomena in the films, and a complete and unifying description of the problem is still not 

■ available. Moreover, some papers devoted to these problem (dynamics of defects) claim contradicting results. These 
contradictions come mainly from the fact that different authors take into account different microscopic dissipation 

"j***; ■ mechanisms, but partially the source of controversy is related to semantics, since different definitions for the forces 
acting on defects are used (see, e.g., discussion in [4]). 

In this paper we examine theoretically the disclination dynamics in free standing hexatic and smectic-C films at 
scales, much larger than their thickness, where the films can be treated as 2d objects. Our investigation is devoted to 
the first (but compulsory) step of defect dynamics studies: a single point disclination in a hexatic or smectic-C film. 
A number of theoretical efforts [5, 6, 7, 8, 9] deal with similar problems. Our justification for adding one more paper 
to the topic is the fact that we did not see in the literature a full investigation of the problem taking into account 
hydrodynamic back-flow effects. Evidently, these effects can drastically modify dynamics of defects. The goal of this 
work is to study the disclination motion in free standing liquid crystalline films on the basis of hydrodynamic equations, 
containing some phenomenological parameters (elasticity modulus, shear and rotational viscosity coefficients). 

In our approach the disclination is assumed to be driven by a large-scale inhomogeneity in the bond or director angle 
providing a velocity for the disclination motion (relative to the film) . As a physical realization of such non- uniform 
angle field, one can have in mind a system of disclinations distributed with a finite density. Then the inhomogeneity 
near a given disclination is produced by fields of other disclinations. One can also think about a pair of disclinations 
of opposite topological charges, then this inhomogeneity is related to mutual orientational distortion fields created 
by each disclination at the point of its counterpart. In fact, the majority of experimental and numerical studies of 
disclination motions in liquid crystals and their models [10, 11, 12, 13, 14, 15, 16, 17, 18] is devoted just to investigation 
of the dynamics of two oppositely charged defects. We solve the hydrodynamic equations and find bond angle and 
flow velocity distributions around the moving disclination. The results enable us to relate the disclination velocity 
and a gradient of the bond angle far from the moving disclination. 

An obvious object, where our results can be applied, is the film dynamics near the Berezinskii-Kosterlitz-Thouless 
phase transition. The static properties of the films near the transition have been investigated in a huge number of 
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papers starting from the famous papers by Berezinskii [19] and Kostcrlitz and Thoulcss [20]. There is some literature 
discussing the theory of dynamic phenomena associated with vortex- like excitations in condensed matter physics: 
vortices in type-II superconductors (see, e.g., [21]), vortices in superfluid A He and 3 He (see, e.g., [22, 23]), dislocations 
in crystals and disclinations (and other topological defects) in liquid crystals (see [10, 11, 12, 13, 14, 24, 25, 26, 27]). 
However, most of the theoretical works on the topic start from phenomenological equations of motion of the defects, 
and our aim is to derive the equations and to check their validity. 

Our paper has the following structure. Section I contains basic hydrodynamic equations for liquid crystalline films 
necessary for our investigation. In Section II we find the bond angle and the flow velocity around the uniformly 
moving disclination in different regions, that allows us to relate the disclination velocity to the angle gradient far from 
the disclination. Different cases, depending on the ratio of the rotational and shear viscosity coefficients, are examined 
in Section III. The last section IV contains summary and discussion. Appendix is devoted to details of calculations 
of the velocity and bond angle fields around the moving disclination. Those readers who are not very interested in 
mathematical derivations can skip this Appendix finding all essential physical results in the main text of the paper. 

I. BASIC RELATIONS FOR HEXATIC AND SMECTIC-C* FILMS 

Let us formulate basic relations needed to describe disclination motion in thin liquid crystalline films. Here we 
investigate freely suspended films on scales, larger than the thickness of the films, where they can be treated as 
two-dimensional objects. Such object can be described in terms of a macroscopic approach, containing some phe- 
nomenological parameters. 

There are different types of the film ordering. Experimentally, hexatic and smectic-C films are observed. In these 
films, as in 3d nematic liquid crystals, the rotational symmetry is spontaneously broken. The smectic-C films are 
characterized by the liquid crystalline director which is tilted with respect to the normal to the film, therefore there is 
a specific direction along the film. In the hexatic films molecules are locally arranged in a triangular lattice, however 
the lattice is not an ideal one. The positional order does not extend over distances larger than a few molecular size. 
Nevertheless, the bond order extends over macroscopic distances. Therefore the phase is characterized by the D^h 
point group symmetry. The order parameter for the smectic-C films is an irreducible tensor Q a p (the subscripts 
denoted by Greek letters run two values, since we treat the films as 2d objects). The hexatic order parameter is an 
irreducible symmetric sixth-rank tensor Q a 0jS^ [28], both such tensors have two independent components. Note that 
the order in the smectic-C films can be readily observed, by looking for in-plane anisotropics in quantities such as 
the dielectric permeability tensor. Because of its intrinsic sixfold rotational symmetry, hexatic orientational order is 
hardly observable. However, it can be detected, e.g., as a sixfold pattern of spots in the in-plane monodomain X-ray 
structure factor, proportional to Qap^s^v (see, e.g., [3] and references herein). 

In accordance with the Goldstone theorem, in the both cases of smectics-C and hexatics the only "soft" degree of 
freedom of the order parameter is relevant on large scales, which is an angle ip (like the phase of the order parameter 
for the superfluid 4 He). It is convenient to express a variation of the order parameter in terms of a variation of the 
angle ip. For the smectic-C films the relation is 

SQafj = -Sipea^Q^p - Sipep^Qa^, (1.1) 

where e ap is the two-dimensional antisymmetric tensor. For the hexatic films 

SQaPySfiu — —5ipe ap Qp0jSnv + • • ■ , (1.2) 

where dots represent a sum of all possible combinations of the same structure. Thus for the both films symmetries the 
order parameter can be characterized by its absolute value |Q| and the phase ip which are traditionally represented 
as a complex quantity '3/ (see, e.g., [24]). For the hexatic films the quantity is written as <J> = \Q\ exp(6z</?), and for 
the smectic-C films the quantity is written as = \Q\ exp(2i(p). 

The angle ip should be included into the set of the macroscopic variables of the films. A convenient starting point 
of the consideration is the energy density (per unit area) pv 2 /2 + e, where p is the 2d mass density, v is the film 
velocity, and e is the internal energy density. The latter is a function of the mass density p, the specific entropy a, 
and the angle ip. In fact, e depends on Vip, since any homogeneous shift of the angle ip does not touch the energy. 
For the hexatic films, leading terms of the energy expansion over gradients of ip are 

e = eo(p,a) + ^(V^) 2 , (1.3) 

where K is the only (due to hexagonal symmetry) orientational elastic module of the film. For smectic-C films 
one should introduce two orientational clastic modules, longitudinal and transversal with respect to the director n, 
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determining a specific direction in the film: Q a p oc 8 a p — 2n a np. However, the anisotropy in real smectic-C films is 
weak because of the small tilt angle of the director. Moreover, fluctuations of the director lead to a renormalization 
of the modules, and on large scales there occurs an isotropization of a smectic-C film [29]. Therefore in this case we 
can use the same expression (1.3) for the elastic energy. 

The complete dynamic equations for the freely suspended liquid crystalline films, valid on scales larger than their 
thickness, can be found in [30]. We are going to treat a quasistationary motion of the disclination. At such motions hard 
degrees of freedom are not excited. By other words, we can accept incompressibility and neglect bending deformations 
which are suppressed by the presence of the surface tension in freely suspended films. In the same manner, for the 
quasistationary disclination motion the thcrmo-diffusive mode is not excited, what implies the isothermal condition. 
For freely suspended films such effects as substrate friction (relevant, say, for Langmuir films) are absent. Thus at 
treating the disclination motion we can consider the system of equations for the velocity v and the angle ip solely. 
The equations are formulated at the conditions p = const, T = const (where T is temperature) and \7v = 0. 

The equation for the velocity follows from the momentum density j = pv conservation law 

d t j a = -V/3 [T aP - v(V a v fJ + Vpva)] , (1.4) 

where T a (j is the reactive (non-dissipative) stress tensor and r\ is the 2d shear viscosity coefficient of the film. For the 
two-dimensional hexatics the reactive stress tensor is (see [30] , chapter 6) 

T a(i = pv ( iv a - qS a fi + KV a ipVpip - y e a7 V 7 V/3<£ - y e^V ~ ( V a (p , (1.5) 

where <r = e — pde/dp is the surface tension. Note that the ratio Kp/rj 2 is a dimensionless parameter which can be 
estimated by substituting 3c? quantities instead of 2c? ones (since all the 2d quantities can be estimated as corresponding 
3d quantities multiplied by the film thickness, and the latter drops from the ratio). For all known liquid crystals the 
ratio is 10~ 3 10~ 4 (see, e.g., [1, 2, 3, 31]), and can be, consequently, treated as a small parameter of the theory. 
The second dynamic equation, the equation for the bond angle, is 

1 K 2 , N 
Ot<p= ~e a pV a vp -v a V a <p-\ V ip, (1.6) 

2 7 

where 7 is the so-called 2d rotational viscosity coefficient. We did not find in the literature values of the coefficient 
for thin liquid crystalline films. For bulk liquid crystals (see, e.g., [1, 2, 3, 31]) the 3d rotational viscosity coefficient is 
usually few times larger than the 3c? shear viscosity coefficient. Therefore one could expect 7 > 77. However, to span 
a wide range of possibilities below we treat r = 7/77 as an arbitrary parameter. 

If disclinations are present in the film, it is no longer possible to define a single-valued continuous bond-angle 
variable ip. However, the order parameter is a well-defined function, which goes to zero at a disclination position. The 
gradient of <p(t,r) is a single- valued function of r which is analytic everywhere except for an isolated point which is 
the position of the disclination. The phase acquires a certain finite increment at each turn around the disclination 

<j> Ar a V a ip = 2ns, (1.7) 

where the integration contour is a closed loop around the disclination position, going anticlockwise, and s is a 
topological charge of the disclination. For hexatic ordering s = (l/6)n, and for nematic orientational order s = (l/2)n, 
where n is an integer number. One can restrict oneself to the disclinations with the unitary charge n (i.e. with s = ±1/6 
for the hexatic films, s — ±1/2 for the smectic-C films) only, since disclinations with larger |s| possess a larger energy 
than the set of unitary disclinations with the same net topological charge, and thus the defects with larger charges 
are unstable with respect to dissociation to unitary ones. Therefore the disclinations with the charges \n\ > 1 do not 
play an essential role in the physics of the films [1, 2, 3, 31]. To bring expressions, given below, into a compact form, 
we keep the notation s for the topological charge. 

The static bond angle is determined by the stationary condition SE/5ip = 0, where 

E = Jd 2 r (V + e) 

is the energy of the film. For the energy density (1.3) the condition is reduced to the Laplace equation 



vV = o. 



(1.8) 
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There is a symmetric solution ip for an isolated motionless disclination which satisfies Eqs. (1.7,1.8), which has the 
following gradient 

V Q ^ - se afj y _ ^ 2 . (1.9) 

Here _R is a position of the disclination. If the origin of the reference system is placed to this point then one can write 
ipo = s arctan(y/x), where x and y are coordinates of the observation point r. If some dynamic processes in the film 
occur, then ip varies in time and the distribution (1.9) is disturbed. It is perturbed also due to the presence of the 
bond angle distortion related to boundaries or other disclinations. 

Below, we have in mind a case, when a system of a large number of disclinations (with an uncompensated topological 
charge) is created. For the 3d nematics it can be done rather easily [1, 2, 3], since due to the intrinsic elastic anisotropy 
the energies of positively and negatively charged defects are different. We do not aware of experimental or theoretical 
studies of defect nucleation mechanisms in free standing smectic or hexatic films. Hopefully, a situation with a finite 
2d density of defects can be realized for the films as well (for instance, the defects could appear even spontaneously 
as a mechanism to relieve frustrations in chiral smectic or hexatic films, i.e. analogously to the formation of the 
Abrikosov vortex lattice in superconductors [32]). Examining a motion of the disclination in this case, we investigate 
a vicinity of the disclination of the order of the inter-disclination distance. Then far from the disclination the bond 
angle ip can be written as a sum const + ur, where u is much larger than the inverse inter-disclination distance (since 
the number of disclinations is large). Near the disclination position the bond angle <p can be approximated by the 
expression (1.9). Our main problem is to establish a general coordinate dependence of ip and v, which, particularly, 
enables one to relate the bond gradient u and a velocity of the disclination. 



II. FLOW AND ANGLE FIELD AROUND A UNIFORMLY MOVING DISCLINATION 

Here we proceed to the main subject of our consideration, which is a single disclination, which is driven by a large- 
scale inhomogeneity in the bond angle ip. A disclination velocity is determined by an interplay of the hydrodynamic 
back-flow and by the own dynamics of the angle ip. To find the velocity, one has to solve the system of equations 
(1.4,1.5,1.6) with the constraint (1.7), requiring a suitable asymptotic behavior. As we explained in the previous 
section, at large distances from the disclination, the angle ip is supposed to behave as const + ur. In addition, we 
are going to work in the reference system, where the film as a whole is at rest. That means, that the flow velocity, 
excited by the disclination, has to tend to zero far from the disclination. 

Let us consider a situation, when the disclination moves with a constant velocity V. In the case both the angle <p 
and the flow velocity are functions of r — Vt (where R = Vt is the disclination position) . Then the equation for the 
velocity (1.4) is written as 



K 

p(Vfj - vp)Vpv a + t]\7 2 v a + —t a (N fi^ 2 <P - K\7 a ip\7 2 ip - 



0. (2.1) 



One can omit the first (inertial) term in the left-hand side of (2.1), which is small because of the smallness of the 
parameter Kp/rj 2 . Then we obtain from Eqs. (1.4-1.5) 

K K 

V 2 «a + ^-ea/jV/jVV Va^VV + V Q w = , (2.2) 

2?7 r] 

where w — 77 _1 [<T — (K/2)(\7p) 2 ]. The equation for the angle ip, following from Eq. (1.6), at the same conditions is 

VV + j(V a V a <p = ^v a V a <p - -Le afj \7 aVl3 . (2.3) 

We are looking for a solution, which is characterized by the following asymptotic behaviour at r — > oo: the velocity 
v vanishes and V</? tends to a constant vector u. It is clear from the symmetry of the problem, that the gradient u of 
the bond angle is directed along the F-axis, provided the velocity is directed along the X-axis. Therefore ip — > uy at 
r — > oo. Our problem is to find a relation between the quantities V and u, that is between the disclination velocity 
and the bond angle gradient far from the disclination. There are two different regions: large distances r 3> u _1 and 
the region near the disclination r <C u^ 1 . At large distances corrections to the leading behavior ip rts uy are small and 
the problem can be treated in the linear over these corrections approximation. In the region near the disclination <p 
is close to the static value (1.9) and the flow velocity v is close to the disclination velocity V (the special case when 
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the ratio 7/77 is extremely small will be discussed in Subsection IIIC). Those two regions are separately examined 
below. The relation between u and V can be found by matching the asymptotics at r <~ u^ 1 . As a result, one obtains 

V=— Cu, (2.4) 
V 

where C is a dimcnsionless factor depending on the dimcnsionlcss ratio r = 7/77. This factor C is of the order of 
unity if r <~ 1. We are interested in the asymptotic behavior of C at small and large r. 



A. Region near the disclination 

Let us consider the region r <C u^ 1 . Here one can write 

f = (p (r - R) + <pi(r - R) , 



(2.5) 



where R = Vt is the disclination position, the angle ipo corresponds to an equilibrium bond angle around an unmoving 
disclination, and ipi is a small correction to ip - The gradients of <po are determined by Eq. (1.9). 
Linearizing the equations (2.2,2.3) in ipi, one obtains 



rf7 2 v 
72 



K 



^-f(v^) 2 



= 0, 



v *y *y *y 

V<Pi - -j^v a V a tp a + —e af:) V a vp = - — V a V a (p a . 

Introducing a new variable \ = (K/i])\7 2 if 1 we rewrite Eqs. (2.6,2.7) as 

2 1 

V v a + -e af3 V(3X ~ Vq^oX + V Q tu = , 

r 

X - rv a V a ip + —e aP W a vp = -rV a d a (p , 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



where as above r — 7/77 and w = 77 1 [<r — K/2(\7Lp) 2 ] . It follows from Eq. (2.8) and \7 a v a — that V 2 tu = 
V Q <poV Q x- A solution of the system (2.8-2.9) can be written as 



v a = V a + tapVpfl , 



(2.10) 



where V a is an obvious (due to Galilean invariance) forced solution and the stream function f2 describes a zero mode of 
the system (2.8-2.9). The system is homogeneous in r, and therefore is a sum of contributions which are power-like 
functions of r. 

Taking curl of Eq. (2.8) we obtain 



- V 4 f> - -V 2 X - e 7tt V Q ^V 7 x = . 



(2.11) 



Substituting \ expressed in terms of v from Eq. (2.9) into Eq. (2.11), and using explicit expressions for derivatives 
of <po, in the polar coordinates (r, <fi) we obtain 



r 



1 



1 



1 



1 + - ] v 4 ^ + s r { —d?n - —v'n - s—dfn + s —d r n 



2 r 



0. 



(2.12) 



Solutions of Eq. (2.12) are superpositions of terms oc r a+1 exp(im^). Substituting this r, 0-dependence into Eq. 
(2.12), one obtains the equation for a, which has the following roots 



a = ± 



V2 



2 + 2m 2 - s(l -s)r±^j (2 + 2m 2 - s(l - s)f) 2 - 4sf (to 2 - 1 + s) - 4(to 2 - 1) 



1/2 



, (2.13) 



where P = T(l + T/4) 1 . Hence < T < 4 for any 7 and 77. Evidently, all the roots (2.13) are real. Let us emphasize 
that there is no solution a = (corresponding to a logarithmic in r behavior of the velocity) among the set (2.13). 
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The first angular harmonic with |m| = 1 is of particular interest because far from the disclination ip\ — ur sin <j> and 
fl = —Vrsmcf). If r is small, there is a pair of small exponents 

a = ±ai, ai = sV7/2, (2.14) 

for m = ±1. Otherwise, for any other relevant m (terms with m — are forbidden because of the symmetry), 
exponents (2.13) have no special smallness. 

We established that fl is a superposition of the terms oc r a+l exp(im^) with a determined by Eq. (2.13). Then the 
velocity can be found from Eq. (2.10). To avoid a singularity in the velocity at small r, one should keep contributions 
with positive a only. By other words, one can say, that the velocity field contains contributions with all powers a 
(2.13) but the factors at the terms with negative a are formed at r ~ a (where a is the disclination core radius), and 
therefore the corresponding contributions to the velocity are negligible at r > o (this statement should be clarified 
and refined for small negative exponents —oti in the limit of small r, see below, Subsection IIIC). We conclude, that 
the correction to V in the flow velocity v, related to f2 in Eq. (2.10), is negligible at r ~ a. Thus we come to the 
non-slipping condition for the disclination motion: the disclination velocity V coincides with the flow velocity v at 
the disclination position. 

Next, to find ip, one should solve the equation (K/rj)V 2 ip = x, where \ is determined from Eq. (2.9). Then, besides 
the part determined by the velocity, zero modes of the Laplacian could enter ipi. The most dangerous zero mode is 
Uy, since it produces a non-zero momentum flux (and associated to it the Magnus force) 

(p <ir a e Q pTpy ~ KU , 

to the disclination core. However, because of the condition a^O, all the contributions to the velocity correspond to 
zero viscous momentum flux to the origin. Consequently, it is impossible to compensate the Magnus force by other 
terms. The above reasoning leads us to the conclusion that the factor U (and therefore the Magnus force) must be 
zero. Thus, <p\ contains only terms, proportional to r a+1 , with a > 0. This conclusion is related to the fact that 
for free standing liquid crystalline films any distortion of the bond angle unavoidably produces hydrodynamic back- 
flow motions (i.e. v ^ 0). Unlike free standing films, for the liquid crystalline films on substrates (Langmuir films) 
hydrodynamic motions (back-flows) are strongly suppressed by the substrate, and a situation when the back-flow is 
irrelevant for the disclination motion could be realized. 



B. Remote Region 

Let us consider the region r 3> u , where we can write ip = uy + (p and linearize the system of equations (2.2,2.3) 
with respect to (p. Then we get the system of linear equations for v and (p: 

K 

V v a + — (e a0 WpV 2 (p - 2u a y 2 <p) + \7 a m = , 

(V 2 + 2pd x )p+ 7 ^(e a pV a vp-2uv y ) = 0, (2.15) 
where p = V-f/2K. Taking curl of the first equation and excluding the Laplacian, we obtain 

epaVpVa = — [(V 2 + 2ud x ) + (2.16) 
where $ is a harmonic function. In terms of the harmonic function $ the system (2.15) is reduced to 

[(i + r/4)v 4 + 2 P v 2 d x - ru 2 dl] <p = {r/2)ud x <$> . (2.17) 

The equation (2.17) can be written as 

(V 2 + 2k l 8 x ){V 2 - 2k 2 d x )$ = {r/2)ud x $ , (2.18) 

kl ' 2 = 2(1 +V/4) (Vp 2 + r(i + r/4)u 2 ± P ) . (2.19) 

The quantities fci and k 2 have the meaning of characteristic wave vectors. We conclude from Eq. (2.18) that zero 
modes of the operator in the left hand side of the equation are proportional to 

exp(— k\r — k\x) , exp(— k 2 r + k 2 x) , 
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that is they are exponentially small everywhere, except for narrow angular regions near the X-axis. The behavior 
of the zero modes inside the regions is power-like in r. Besides, there is a contribution to Cp related to the harmonic 
function <!>. It contains a part which decays as a power of r (the leading term is oc r _1 ) at large r » u^ 1 . This 
solution is examined in more detail in Appendix A 1 . 

III. DIFFERENT REGIMES GOVERNED BY r 

A character of the velocity and the bond fields around the moving disclination is sensitive to the ratio of the 
rotational and the shear viscosity coefficients r = 7/77. In this section we examine different cases, depending on the 
r value. 

A. The case r > 1 

We start to treat different regimes of mobility with the most likely case r > 1. If r ~ 1 then the factor C in 
Eq. (2.4) is of order unity and u <~ p. Then we find from Eqs. (2.19), that fci, fc 2 ~ u. It is a manifestation of the 
fact, that there is the unique characteristic scale in this case, which is u . Then one can estimate (p by matching at 
r ~ of the solutions in the regions near the disclination and far from it. We conclude, that it is a function of a 
dimcnsionless parameter ur; the function is of the order of unity, when its argument ur is of the order of unity. 

For large r there remains the unique characteristic scale w" 1 and, consequently, in this case still C <~ 1. To prove 
the statement, we treat first small distances r <C u~ x . As we have shown in Subsection II A, the corrections <p\ and 
Sv, to (po and V respectively, are expanded into the series over the zero modes characterized by the exponents (2.13). 
Particularly, for m = 1 we can write ipi ~ uy(ur) a . In the large r limit the exponents a (2.13) are regular since 
r — > 4. From (2.13) we have ct\ ~ 1, and for this case 

K 

X ^uy(ur) ai . 

Comparing equations (2.8) and (2.9), we conclude that at large F the term with \ m tne second equation can be 
omitted and therefore the equation is a constraint imposed on the velocity. Then the first equation gives 

K 

\5v\ ~ — uy(ur) ai . 
r\r 

The disclination velocity can be found now from the relation V <~ \Sv\ at the scale u -1 , that is p <~ Ax, or C ~ 1. 
The complete analysis covers also the remote region. At the condition p ~ _Tm both ki, 2 ~ u^ 1 . Then using the 
procedure of Appendix Al, one can prove that the solutions from two regions can be matched at r ~ u^ 1 , and 
thus there are no new characteristic scales, indeed. We also note that the rotational viscosity 7 is excluded from 
the hydrodynamic equations at large _T. Although this is not true inside the disclination core provided that one can 
employ hydrodynamics there (see Appendix A 4), boundary conditions for v and <p on the core reveal no dramatic 
behavior. Consequently it is the shear viscosity alone that determines the disclination mobility, which means that 
C~ 1. 

Thus we can say that the limit r — > 00 is trivial: no any additional feature appears in comparison with r ~ 1. 
However, it is not the case for small r, since for r <C 1 it turns out that u 3> p. We study this case in the next 
Subsection. 



B. Small r 

Here, we consider the case J 1 « 1. This limit is attained at anomalously large 77, so that Kp/rj 2 should be still 
treated as the smallest dimensionless parameter. That justifies use of the same equations (2.2,2.3), as in the previous 
subsections. 

At r <C u^ 1 the analysis of Subsection II A is correct. As we noted, contributions to v and ipi, related to the modes 
with negative a, should not be taken into account there. At r <§C 1 the leading role is played by the mode with the 
smallest exponent (cti = s\/T'/2), because the presence of modes with positive exponents a <~ 1 would contradict the 
condition of smooth matching at r ~ u" 1 . Strictly speaking, neglecting a small negative exponent —ot\ is correct at 
the condition a\ \ ln(ua)| » 1, where a is the core radius of the disclination. In this subsection we consider just this 
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case. The opposite case, that we will term it as extremely small r limit, is analyzed in Subsection IIIC. At r <C u 1 
we can write 



if! ~ uy(ur) ai , 



V — v x <~ aiu— (ur) Ql , 

7 



(3.1) 



where the coefficient at y(ur) ai is determined from matching at 



where V93 ~ 1/r. Analogous matching 



U — f 1 ~ ^ at r ~ m 1 gives V ~ aiuK/j. The relation can be rewritten as p ~ aiu <C w. Therefore we conclude 

c~ i/Vr. 

The relation p ~ vT' u leads to k\ , 2 ~p<u, as follows from Eq. (2.19). By other words, a new scale (different 
from m _1 ) appears in the problem. Therefore a detailed investigation of the far region r> u' 1 is needed, to establish 
an r-dependence of the bond angle tp and the velocity field v there. This investigation can be based on the equations, 
formulated in Subsection II B, which are correct irrespective to the value of pr. 

Explicit expressions, describing the velocity and the angle, are presented in Appendix Al. They contain three 
dimensionless functions £i(V/u), c i(V/w) and c 2 (V/u). Really, at ur ^> 1 one can keep only zero terms of the 
expansions of these functions into the Taylor series. Among these three coefficients only one is independent, see Eq. 
(A10). Therefore the general solution could be expressed with a single free parameter, ( = (i(0). The procedure 
corresponds to the following construction of solutions of the equations of motion (2.15) in the region ur S> 1. We 
have to match the solutions in the outer (far from the disclination) and in the inner (close to the disclination) regions 
at ur ~ 1. Technically the matching is equivalent to the appropriate boundary conditions for the outer problem at 
ur ~ 1, and these boundary conditions can be replaced formally by the local source terms in the equations, acting 
at ur ~ 1. One can expand these sources in the standard multipolar series. Thus we come to the expansion over 
the gradients of S- function. The gradients scale as u, and therefore £, Ci, c 2 are dimensionless functions of the 
dimensionless ratio V/u. 

To establish asymptotic behavior of the angle ip and of the velocity v, let us first consider the region u^ 1 «r« 
Then we derive from Eqs. (A4,A5,A10) 



K(2s-C) 
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k\k 2 ki(pr) . 



Here we presented only the leading logarithmic contribution of the zero harmonic in v x . 



derivatives, determined by Eqs. (3.1,3.2), at r 
(A2,A5,A10), we obtain for the region u 



1 < r ^p- 1 



we find C ~ 1 (we imply s ~ 1). 



ip = ipa + uy + spy ln(pr) . 



(3.2) 

Matching the velocity 
Using the expressions 

(3.3) 



We see, that there is only a small correction to the simple expression (p + uy there, since p^u. 

In the region pr >• 1, the expressions for the angle ip and the velocity v are more complicated. Using the formulas 
(A2,A3,A4,A5), one obtains 



dxP - 



ITT 

Hj 2 



C\\fk[ exp(-fcir - k\x) + ci\fk~2 exp(-fc 2 r + k 2 x) 



dyip = u — 2s J — c\\fk\ cxp(— k\r — k\x) — ci\fk~2 cxp(— k 2 r + k 2 x) 



y _ (y_ 

r 3/2 2 r 2 ' 
1 C x 

H72 + 2^2' 



"fU 



c\k 2 \fk~\ exp(-fcir - fcix) - c 2 k\\fk~ 2 exp(-fc 2 r + k 2 x) 



y 

r 3/2 



7 u \ V 2 



. exp(— k±r — k\x) H — exp(— k 2 r + k 2 x) 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



where C\ ~ 1 and c 2 ~ 1 are determined by Eq. (A10) (we omitted the argument to simplify the notations). The 
expressions (3.4,3.5,3.6,3.7) contain two types of terms: isotropic and anisotropic ones. The anisotropic contributions 
are essential only in the narrow angular regions near the X-axis, where they dominate. It is worth to note the very 
non-trivial structure of the flow, where the isotropic flux to the origin is compensated due to the anisotropic terms. 

The expressions found in this subsection generalize the famous Lamb solution for the hydrodynamic flow around a 
hard cylinder, (see, e.g., [33, 34, 35]) where far from the cylinder the velocity field is exponentially small everywhere 
except in the wake of the corps, i.e. in the very narrow angular sector ( "tail" ). Disclination motion in liquid crystalline 
films can be regarded as motion of a cylinder framed by a "soft" (i.e. deformable) orientational field ip. Due to the 
additional (with respect to the classical Lamb problem) degree of freedom, our solution has two tails around the 
moving disclination: wake, beyond the disclination, and precursor in front of it. In fact, both degrees of freedom (the 
flow velocity and the bond angle) are relevant. 
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C. Extremely small r 

At the above analysis we implied the condition a\ \ ln(ua)| » 1 (remind that a\ — at small r), imposing a 

restriction from below to r at a given u. If oti\ ln(ua)| <C 1, then the terms with both a = ±ai, determined by Eq. 
(2.14), should be taken into account near the disclination, that leads to a logarithmic behavior of the correction ipi 
to (po there: 

ip i <~ uy\n(r/a)\ ln(aw)| _1 , (3.8) 

instead of Eq. (3.1). Matching the derivatives of expressions (3.3,3.8) at r ~ w _1 gives p ~ u|ln(cra)| _1 . By other 
words, C ~ [-T ln(au)] . This case corresponds formally to the limit r\ — > oo in our equations, when one can drop 
the back-flow hydrodynamic velocity in the equation for the bond angle. The situation was examined in the works 
[6, 7, 8, 9]. We present a simple analysis of the case in Appendix A 2. Note also that, really, there is no crossover at 
r ~ u^ 1 in the bond angle behavior in this situation. 

Let us clarify the question, concerning the Magnus force in this case. In accordance with Eq. (3.8) the reactive 
momentum flux to the disclination core is 



<j> dr a e a pTfj 7 <~ Ku ln(r/a) | ln(au) 



Thus, the flux is r-dcpcndcnt, tending to zero at r — > a. This reactive momentum flux is compensated by the viscous 
momentum flux, (related to derivatives the flow velocity v), which is non-zero in this case because of the logarithmic 
in r behaviour of the flow velocity near the disclination. The flow velocity can be found from Eqs. (2.6,3.8), it behaves 
as 

Va m — \~\ W V /3 b ln {r a)} , 

77|m(ait)| 

what is a generalization of the Stokes-Lamb solution [33, 34]. However unlike the Lamb problem (a hard cylinder 
moving in a viscous liquid), in our case \V — v(r = a)\ <~ V, i.e. we have a slipping on the core of the moving 
disclination. This slipping seems natural in the limit of extremely small values of r, corresponding to the limit 
7] — > oo, that is to a strongly suppressed hydrodynamic flow. Physically, the property means that the disclination can 
not be understood as a hard impenetrable object in the case. It is worth noting also that found above logarithmic 
behavior looks similar as a general feature of the two-dimensional hydrodynamic motion, which comes from the well 
known fact (see e.g., [33, 34, 35]) that for two-dimensional laminar flow even for a small Reynolds number one cannot 
neglect nonlinear terms, which become relevant for large enough distances. However, in our case these non-linear 
terms do not come from the convective hydrodynamic nonlinearity, they come from the nonlinear over <p terms in the 
stress tensor (1.4). 

An explicit expression for tp and its asymptotics, corresponding to the considered case, are found in Appendix A 2. 
An expression for the flow velocity field, excited by the disclination motion, is derived in Appendix A3. 



IV. CONCLUSION 



Let us sum up the results of our paper. To understand physics, underlying the freely suspended film dynamics, 
we studied the ground case, namely, a single disclination motion in a thin hexatic (or smectic-C) film, driven by 
a inhomogeneity in the bond (or director) angle. We investigated the uniform motion (i.e. one with a constant 
velocity). For this case we derived and solved equations of motion and found bond angle and hydrodynamic velocity 
distributions around the disclination. It allows us to relate the velocity of the disclination V to the bond angle gradient 
u = |V<p| in the region far from the disclination. It is in fact the reason, why so much efforts are needed: the full set 
of the equations should be solved everywhere, not locally only. We established the proportionality coefficient C (see 
Eq. (2.4)) in this non-local relationship which has a meaning of the effective mobility coefficient. The coefficient C 
depends on the dimensionless ratio r of rotational 7 and shear viscosity 77 coefficients. 

There is little experimental knowledge as far as the values of the coefficients 7 and r\ in hexatic and smectic-C films 
are concerned. It is believed generally that the corresponding values in a film (normalized by its thickness) and in a 
bulk material are not very different [3, 31], and in this case we are in the regime of F ~ 1, when the coefficient C 
is of the order of unity. However f<l case is not excluded from the both theoretical and material science points 
of view. We found for small fcl case the coefficient C ~ 1/y/T. We established a highly non-trivial behaviour of 
the flow velocity and of the bond angle, which is power-like in r near the disclination and extremely anisotropic far 
from it. And only for the case of extremely small r< 1/ ln 2 (ua) (where a is the disclination core radius), we found a 
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logarithmic behavior C ~ \T ln('ua)] . The main message of our study is that the hydrodynamic motion (that is the 
back- flow), unavoidably accompanying any defect motion in liquid crystals, plays a significant role in the disclination 
mobility. Experimental evidence (see, e.g., the recent publication [36]) shows that it is indeed the case. 

Our consideration can be applied to the motion of a disclination pair with opposite signs. In this case the role by 
the scale u -1 is played by the distance R between the disclinations. Then we find in accordance with Eq. (2.4) that 
dtR oc R" 1 without a logarithm (provided the twist viscosity coefficient 7 is not anomalously small, see Subsection 
IIIC for the quantitative criterion). This conclusion is confirmed by results of numerical simulations for 2d nematics 
[15, 16, 17, 18]. The authors of the works consider the equations of motion in terms of the tensor order parameter, 
consistently taking into account the coupling between the disclination motion and the hydrodynamic flow. They 
computed dynamics of the disclination pair annihilation and found that the distance R between the disclinations 
scales with time t as i 1 / 2 , without logarithmic corrections (as it follows from our theoretical consideration for any 
(excepting extremely small) value of the parameter r. Unfortunately we did not find in the publications [16, 17, 18] 
the magnitudes of the the shear viscosity which was used in the simulations. Lacking sufficient data on the values of 
7 and r\ we can at present discuss only the general features of the disclination dynamics. For instance the authors 
[18] found numerically in one-constant approximation the asymmetry of the disclination dynamics with respect to the 
sign of the topological charge (s = ±1/2). In our approach the asymmetry naturally appears from nonlinear terms in 
the stress tensor (1.5) and from the first term in the right-hand side of Eq. (1.6) responsible for the different coupling 
of orientational and hydrodynamic flow patterns for positive and negative disclinations. This results in the fact that 
smaller positive exponents in Eq. (2.13) (corresponding to the minus in the brackets) for each m are larger if s = 1/2 
than those for s = —1/2. Thus, the disclination with s = 1/2 exerts a stronger influence on the flow velocity, the 
qualitative conclusion that was arrived in [18]. 

Although the theory, presented in the paper, is valid for free standing hexatic or smectic-C liquid crystalline films, 
the general scheme can be applied to Langmuir films, that is for the liquid crystalline films on solid or liquid substrates. 
Since the Langmuir film is arranged on the substrate surface, then any its hydrodynamic motion is accompanied by 
the substrate motion. Therefore for solid substrates a situation, when the hydrodynamic back-flow is irrelevant for 
the disclination dynamics, could be realistic. In Subsection IIIC (see also Appendix A 2) we examine this limit, 
and reproduce the results of the works [6, 7, 8, 9], where the hydrodynamic back-flow was neglected from the very 
beginning. As it concerns to the Langmuir films on a liquid substrate, this case requires a special investigation, 
however the approach and main ideas of our paper could be useful as well. 

Our results can be tested directly by comparing with experimental data for smectic-C films. The hexatic order 
parameter, which has sixfold local symmetry is not coupled to the light in any simple way (and therefore ideal hexatic 
disclinations are hardly observed in optics). However it is possible to observe the core splitting of the disclinations in 
tilted hexatic smectic films [26] . Indeed, due to discontinuity the tilt direction (which is locked to the bond direction) 
one can observe hexatic order and hexatic disclinations indirectly. The second possibility of checking our theoretical 
results is the classical light-scattering (where in typical experiments wave vectors are q = 10 2 -j- 10 cm , and the 
frequency is u> < 10 8 s _1 ). For a reasonably thick film the power spectrum of light scattering can have some additional 
structure that reveals the disclination properties (e.g. defects are thought to be relevant to the very low-frequency 
noise observed in thin films). Such kind of experimental studies are highly desirable. 
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APPENDIX A 

Here we present some technical details, concerning the velocity and the bond angle fields around the disclination. 
The case of small r is mostly treated, which is especially rich from the point of view of the space structure of the 
fields. 
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1. Distances far from the disclination 



Here we derive some results for the far from the disclination region which are afterwards applied to the case of 
small r considered in Subsection III B. 

Let us examine the harmonic function $ in Eqs. (2.16). Since the function is analytic in the region r > it can 
be expanded over derivatives of lnr there. Next, due to the symmetry of the problem, $ should be antisymmetric in 
y. Therefore at least one derivative d y should be present in each term of the expansion, that is 



$ = u(id y lnr , 



(Al) 



where (i = Ci(V/u), and (i(z) is a series of z } converging in a circle with radius of the order of unity. The expansion 



coefficients in the series (i (V / u) are determined by matching at r 



with the inner problem. 



Because of the symmetry, the angle <p can be represented in the following form: 

d x (p = d y B, 8 y ip = -{H + d x B), \7 2 B + d x H = 0. 



(A2) 



The latter equation is the condition t a pSJ ' Q V p<p = 0. Note that V 2 (p = —d y H. In the far from the disclination region 
we can use Eqs. (2.15,2.16). The incompressibility condition V Q u Q = must be taken into account as well. Then one 
obtains expressions for the velocity in terms of B and H: 



curl v = — d v 
2r\ v 



-H + 2uB + u(i ln(pr) 



v x = 



K 

ju 
K 
-fu 



H + 2pB 



r r 



d x { -H + 2pB + 



r 



H + 2uB + u(i ]n(pr) 
-H + 2uB + u(x ln(pr) 



~2~ 



-H + 2uB + u(i ln(pr) 



(A3) 
(A4) 



Solutions of Eq. (2.17) imply: 



B 



^Koih^e-^ + c 2 K (k 2 r)e 



k 2 x 



2& Mpr) ■ 



H = 2s [k^Koik^e-^ - k 2 c 2 K (k 2 r)e 



k 2 xl 



(A5) 



Here the particular representation (Al) is used and an arbitrary function of y which can contribute to H is chosen 
to be zero because Vt^ — > (and hence H — > 0) as r — > oo. Above t\ and c 2 are dimensionless differential operators 
which can be represented as Taylor series over V/it, i.e. ci(V/u) and c 2 (V/u). These functions must scale with u 
since the functions have to be found from the matching at r ~ u^ 1 . 

Additionally, there are two conditions on the variables in the region tir > 1. First, the correct circulation around 
the origin leads to the effective 5- functional term in Eq. (A2): 



X7 2 B + d x H = -2ns6{r). 
The second condition is the absence of the flux to the origin: 

d(j)v r (r, <j>) = . 



The relations (A6,A7) lead to the conditions 
ci(0) + c 2 (0) + Ci(0)/2s = 1, 



/• 



i 



[fcici(0)-fc 2 c 2 (0)] - [p + 



ru 

[d(0)+C2(0) + Ci(0)/2*] + — Ci(0) 



At small r the solution of the equations (A8,A9) is 

Ci(o) = c, ci(o) 



ki - Ck 2 /2s 



k\ + k 2 

We also assumed that Q < 1, which is justified in Subsection IIIB. 



c 2 (0) = 



k 2 - Cfci/25 
h + k 2 



(A6) 

(A7) 

(A8) 
(A9) 

(A10) 
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2. Suppressed Flow 

Here we demonstrate, for convenience of references, how is it possible to find the disclination velocity V provided 
the hydrodynamic velocity v is negligible (say, due to a substrate friction). We reproduce the results of the works 
[6, 7, 8, 9]. 

In the absence of the hydrodynamic flow the equation for the angle p is purely diffusive: 

jd t p = KV 2 p, (All) 

as follows from Eq. (1.6) at v = 0. We assume that <p — > uy if r — > oo. The disclination motion is forced by the 
"external field" u. Next, we are looking for a solution <p(t, x, y) — <p(x — Vt, y). Then we get from Eq. (All) 

2pd x <p + V 2 <p = 0, where 2p = -)V/K. (A12) 
Below we consider a solution corresponding to a single disclination with the circulation 

drV(p = 2-n-s, (A13) 



where the integral is taken along a contour, surrounding the disclination in the anticlockwise direction. The quantity 
s in Eq. (A13) is an arbitrary parameter (which is equal to ±1/6 for hexatics and ±1/2 for nematics). For a suitable 
solution of Eq. (A12), corresponding to Eq. (A13), we have 

/d^ 1 
exp(iqr) = s exp(-px)d y K {pr) . (A14) 
2ir q A - 2ipq x 

This derivative tends to zero at r — > oo as it should be. 

The expression (A14) does not determine tp itself unambiguously, since d x (uy) — and therefore one can add an 
arbitrary function uy to a solution and it gives a new solution then. Note that uy is a zero mode of the equation 
(A12). Thus the solution can be written as 



/•OO 

ip = <f) + uy, <t>{x',y) = -s / dx cxp(-px)d y K (pr) , 

J X 1 



(A15) 



where <f> tends to zero if r — > oo. To relate u in Eq. (A15) and p, one must know boundary conditions at r — > 0, really, 
at r ~ a, where a is the core radius. At small r (p can be written as a series (p = ipo + tpi + . . . , where po corresponds 
to an motionless disclination, p 1 is the first correction to p , related to its motion. Matching with the inner problem 
gives 

Vpi(a)~p, (A16) 

since the solution for the order parameter inside the core is an analytical function of r/a, and the expansion over p is 
a regular expansion over pa. 

Expanding Eq. (A14) over p we get at pr <C 1 

1q y , vxy 

-d x <p W ^ + "a" • 

Therefore in accordance with Eq. (A15) one gets with the logarithmic accuracy 

pi = spyhi{pr) + uy . (A17) 
Using now the boundary condition (A16) one gets with the same logarithmic accuracy 

u = spin (—] . (A18) 



pa 



It can be rewritten as 
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The same answer (A19) can be found from the energy dissipation balance. First of all, we can find the energy E 
corresponding to the solution (A15): 



E = J d 2 ry(Vrf^ J d 2 r 



l„2 + I(v^) 2 + «a B 



(A20) 



Here the first term is the energy of the external field and the second term is the energy of the disclination itself. 
Obviously, both the terms are independent of time, and only the last cross term depends on the time. For \x — Vt\ » 



P- 1 



if x > Vt , 
-2tts if x < Vt . 



/oo 
dyd y (f>^ 
-oo 

Therefore we obtain from Eq. (A20) 

d t E = -2swKuV, (A21) 
On the other hand, we can use the equation (All) to obtain 

dtE = -— ( d 2 r (VV) 2 • (A22) 



7 

Substituting here V 2 ip by 2pd x tp in accordance with Eq. (A12) we get 

OtE^-jV 2 J d 2 r(c^) 2 . 

The main logarithmic contribution to the integral comes from the region a < r < p^ 1 where d x ip w ~sy/r 2 . Thus we 
obtain 

dtE = -irs 2 jV 2 In . (A23) 

Comparing the expression with Eq. (A21) we find the same answer (A19). 

3. Extremely small r 

Here we consider the flow velocity, excited by the moving disclination in the case of the extremely small r. The 
velocity is zero in the zero in r approximation (the case is considered in Appendix A 2), therefore we examine the 
next, first, approximation in r. We use the same formalism and the same designations as in Appendix A 1. 

In accordance with Appendix A 1 solutions of the complete set of nonlinear stationary equations can be represented 
in the following form: 

d x (p = d y B, d y (p= -(H + d x B), (A24) 
K 

curl v = — [-d y H + 2ud v B + 2usd y In r + $'] , (A25) 
lr\ 

v x = -^-dyX/- 2 [-d y H + 2ud y B + 2usd y In r + $'] , (A26) 
2rj 

v v = ^ X V -2 [~d y H + 2ud y B + 2usd y \nr + $'] , (A27) 
Zrj 

where B, H and $' are to be found from the equations 

- d y H + 2pd y B + jV" 2 (V 2 - 2ud x ) [-d y H + 2ud v B + 2usd y Inr + = 

- ^ (dyX/- 2 [-d y H + 2ud y B + 2usd y Inr + $'] d y B + ^V" 2 [-d y H + 2ud y B + 2usd y Inr + $'] (d x B + H^j , 

(A28) 
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*' - 2\7- 2 [(d x B + H)d x d y H + d y Bd 2 y H] , (A29) 
V 2 B + d x H = -2nsS(r) . (A30) 

If r is extremely small, s 2 r\v?(ua) <C 1, the solution of Eqs. (A28-A30) can be continued to the vicinity of the 
core. In the main approximation the solution for ip coincides with the solution for the angle (f>L in the absence of 
the liquid motion. This case, examined in [6, 7, 8, 9], is described in Appendix A 2. The functions Bl and Hl, 
corresponding to tfiL, are given by 



2pB L =H L = 2spK (pr) exp(-px) . (A31) 



This solution results in 



* = 2^»ln(=*^). (A32) 
Then, after neglecting the nonlinear right-hand side of Eq. (A28), we can find 

H[T) ~ 1 + r/4 J (2nf CMlQr) (q 2 2ik iqx )(q 2 + 2ik 2 q x ) ' (A33) 

B{r) can be found in the similar way. Using B and H, from Eqs. (A26,A27) we calculate such flow velocity v(r) that 
vanishes at infinity. 

For r >• p^ 1 this solution coincides with expressions (A5,A8,A9) with 

Ci(0) = 2 S + ^lnf- 
u \pa 

For pr < 1 the expression (A33) is reduced to (A31) and this region produces the main contribution to $' (A32). 
The following expressions are obtained in the inner region (pr <ti 1) from the solution (A24-A33): 

w\= \u-sph\—)y + spy ln - , (A34) 
\ pa J a 

cm \ v ^^l\ n L v. (A35) 
rj a r z 

A relation between p and u is fixed by the condition (A16), leading to u = spln[l/(pa)], which is equivalent to Eq. 
(A19). The flow velocity at pr <§C 1 and ln(r/a) 3> 1 is 

s 2 r 

v a = —Ve a0 Vp [y In 2 (r/a)] , (A36) 

which corresponds to the stream function 

n _-V,-^M.®. (A3T, 



4?7 

The expansion over r near the disclination is regular and can be derived from Eqs. (2.8,2.9) with the condition 
V(fi (a) <~ p: <pl + uy is the zero term of the series for ip, and expression (A37) represents the zero and the first terms 
for 0. 

Note that in accordance with Eq. (A36) in the limit P — > the flow velocity tends to zero near the disclination 
core, v(a)/V = 0(r), despite the fact that the disclination itself moves with the finite velocity V, thus there is a 
slipping on the disclination core in this limit. 



4. Solution with the complete order parameter 



Here we consider the dynamic equations for the coupled velocity field v and the complete order parameter ^ = 
Q exp(iip/\s\) of the hexatic or smectic-C films. These equations are needed to examine the velocity field in a vicinity 
of the disclination core. Below, we work in the framework of the mean field theory. 
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Formally the equations can be derived, using the Poisson brackets method [30, 37]. The energy, associated with 
the order parameter, in the mean field approximation is 



"•4/ A ( |B| '^ ,i - |,|,) 1' 



its density passes to the Af-contribution in Eq. (1.3) on large scales r » a s . The only non-trivial Poisson bracket, 
which has to be added to the standard expressions (see [30, 37]), is [28] 

O'a(ri), *(r 2 )} = -V a **(ri - r 2 ) + ^-r *(r 2 )e Q ^V /3 ^(ri - r 2 ) . 

Z \S\ 

The expression for the energy and the Poisson bracket is correct for hexatics, and turns to the above expressions at 
large scales for the smectics-C (due to fluctuations). Then the dynamic equations read: 

pd t v a + pv fj V fj v a = V V 2 v a - jv a f* (V* + ^* (1 - l*| 2 )) + (v 2 ** + ^* (1 - |lf)) | 

- ^c af 3Vp {**V 2 * - *V 2 **} + V a S, 
i T „ ^s 2 /„9 T 1 



9** + v a V a * = -^*e Q/3 V aV/3 + ^- ( V 2 * + — * (1 - |*| 2 ) , (A38) 
Z\s\ Z--f s \ a s j 

the relation 7 S = s 2 -f/2 ensures the reduction to Eq. (1.6) in the limit \ty\ = 1, and the kinetic coefficients are believed 
to be independent of Q (otherwise one can assume, for example, the dependence 7 S = s 2 -f\^\ 2 /2). The slow dynamics 
of a 2d liquid crystalline system with disclinations can be described by Eqs. (A38) with the additional condition of 
incompressibility Vi> = which enables to exclude the passive variable <f. 

If the distance from the disclination point to another disclination is much larger than a s , i.e. the perturbation to 
the static solution \?o = Qo exp(i<po/|s|) for the single defect is small, we can linearize Eqs. (A38) with respect to the 
perturbation expressed in terms of corrections Q\ and ipi to Qo and tpo respectively: 

r]\7 2 v a - 2 7s jv Q Qo(v/3 - Vp)VpQ + ^QgV a po ((vp - Vp)Vp<po - \e 01 V pv^j J (A39) 



V Q c - , 



^ (v 2 Ql - ^ff Ql - 2V ^ 2 V ^ Q + \ (1 - 3Q 2 ) Ql ) = («„ - Vp)Vp Qo , (A40) 

Ks 2 ( „ \ 1 

V 2 ^i + 20o 1 (V ce QiV Q ^o + V a Q V Q ^i) ) = --e aP V a v p + (v p -Vp)Vp<po. (A41) 



2 7s V v^-v.-^r. ~-^.<*t.,j 2 

In terms of dimcnsionless quantities L = r]Q/K, R = r/a s and r = 2'-f s /(s 2 r/), Eq. (A39) takes the following form 
(as previously, we consider a disclination with the unitary topological charge |s| or — |s|): 

<L + r r 2 QBQ£_p + ( v2r+ 2s \ / Q , / ^ _ v\ A) = 0; (A42) 



<5o is found from 

d 2 R + ^d R -^jQ + Qo(l-Ql) = 0, Qo(0) = 0, Q (oo) = l. 

If i" 1 3> 1, as it follows from Eq. (A42), a new scale R ~ 1/vT <C 1 appears inside the core, the first term in Eq. 
(A42) can be neglected on larger scales and there is no crossover at R ~ 1. 

If Qo = 1, Eq. (A42) is reduced to Eq. (2.12). If R < 1, Q = AR (A w 0.58), and Eq. (A42) can be rewritten as 

V 2 Jv 2 f L+^(i? 2 V 2 ? -4 S 2 )4=0. 
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The solution of the equation is a composition of terms A(i?) sin(m<^>) with different m. The condition A(i?) = keeps 
two constants in the general solution of the ordinary differential equation for X(R), two regular near R = partial 
solutions are 

*M and *M a * f M - HWgW 1 + > 



where 2-F1 is the hyper-geometric function ( 2 Fi(a, 6, c, z) = 1 + a6z/c+ . . . ). Two constants (e.g. derivatives A^' m ^(0) 
and A(' m ' +2 )(0)) are chosen to ensure the slowest possible grow at R 3> 1, i.e to exclude the largest exponent among 
a in Eq. (2.13). 

If r » 1, it is possible to derive a better approximation in the core. We can expand Qo(R) in a series, look for 
the series solution \{R) and extract the terms of the highest order in r. For example, for m = 1 the series for \(R) 
begins with liR + I3R 3 , which fix two constants in the partial solution: 



X(R) =hR 
+h 



1 + 



A 2 r s (2 - s 2 ) 



s 2 A 2 TR 2 1 1 - VI + 4s 2 1 + VI + 4s 2 
o 1" 2-F1 7, > 7. 



, 2, 



A 2 TR 2 



A 2 Ts 2 



R 



-1+ 2^1 



1 - Vl + 4s 2 1 + Vl + 4s 2 



, 2, - 



A 2 ri? 2 



The solutions of Eqs. (A40,A41) have the following form: 

Qi — , &{R)d ( f l sin(m</>) , ipi — a(R) sin(m</>) , 
and should be found from the equations 



1 



i? 2 ' 



i? 2 



■1? + d R Q <j' 



r 

~2 



A" + 



which generalize the expressions of Ref. [7] . 

The dynamic equations with the complex order parameter demonstrate that at any r nothing special happens 
on the core in the boundary conditions for Eqs. (1.4-1.6). The peculiarity of extremely small r which leads to 
non-slipping condition is in the slow grow of Vfi far from the disclination. 
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